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Abstract. We investigate the phonon thermal conductivity « p h of doped La2Cu04 based on out-of-plane 
thermal conductivity measurements. When room temperature is approached the temperature dependence 
of K p h strongly deviates from the T _ """-decrease which is usually expected for heat transport by acoustic 
phonons. Instead, K p h decreases much weaker or even increases with rising temperature. Simple arguments 
suggest that such unusual temperature dependencies of K ph are caused by heat transport via dispersive 
optical phonons. 

PACS. 66.70.+f Nonelectronic thermal conduction and heat-pulse propagation in solids; thermal waves - 
74.72.Dn La-based cuprates - 44.10.+i Heat conduction 



The thermal conductivity n is an interesting tool in 
order to probe dissipation and scattering of any propa- 
gating excitation in a solid. A recent example is k of com- 
plex materials with low-dimensional spin structures where 
magnetic excitations provide an unusual transport chan- 
nel of heat PEEEBS]. Such magnetic heat conduction can 
in general only be measured in parallel with the phononic 
heat transport of the underlying crystal lattice. A thor- 
ough understanding of the phonon thermal conductivity 
K p h is therefore essential to identify and separate a mag- 
netic contribution K mag - 

In early experiments on magnetic materials clear de- 
viations from a low-temperature scaling as T 3 , which as 
T — > is expected for K p h, were assessed as one im- 
portant indication of K mag [SEIHl- Recently, experiments 
have been performed on materials where magnetic cou- 
pling and velocities of magnetic excitations are several 
orders of magnitude larger than in these early studies. 
Significant K mag is in these cases present at much higher 
temperatures where K p h is expected to follow a T _1 be- 
havior. Prominent examples are the spin-ladder system 
(Sr, Ca, La)i4Cu2404i and the two-dimensional antiferro- 
magnetic La2CuC>4, where pronounced peaks are observed 
in k at high T [111211511115] . These high-T anomalies reflect 
the dimensionality of the underlying magnetic structure, 
i.e., they are only observed when k is measured along the 
direction of ladders or, in the case of La2CuC>4, parallel 
to the magnetic planes. However, it is not a priori clear 
that such anomalies result from magnetic heat transport. 
Just as well, it is possible that dispersive optical modes 
provide a thermal transport channel which could gener- 



ate an unusual high-T behavior of K p h. Even the observed 
anisotropics could be explained in such a scenario, since 
low-dimensional magnetic structures often originate from 
lattice sub-structures with a similar low dimensionality. 
It is therefore essential to carefully investigate the origin 
of high T anomalies. A firm proof that indeed magnetic 
heat conduction is present, is for example obtained if the 
analysis of the additional contributions yields characteris- 
tic magnetic properties like the spin gap of the ladders or 
magnetic defect distances Q]2|5J. 

The actual relevance of heat conduction by optical 
phonons has, however, never been studied systematically. 
It is the purpose of this paper to initiate future work 
in this field by pointing out that optical phonons pro- 
vide a substantial transport mechanism of heat in doped 
La2CuC>4. Our estimation for the heat conductivity by 
optical phonons K op t is based on realistic values for the 
phonon energies and velocities in this material. We find 
that at room temperature re op t can amount to about 40% 
of the total phonon thermal conductivity K p h- 

The experimental data which give reason for our con- 
siderations is the out-of-plane thermal conductivity k c of 
Lai.8-xEuo.2Sr x Cu04-single crystals with x = 0, 0.08, 
0.15 and 0.2. Fig. ^ displays k c as a function of T. Prior 
to discussing the high-T behavior of k c in more detail, 
we briefly summarize the most important facets of the 
doping- and temperature dependence of k c which have 
been elaborately discussed in Ref. Details of the ex- 
periment are also described there. Unlike the in-plane- 
thcrmal conductivity of doped La2Cu04 which, de- 
pending on doping, consists of phononic, magnetic and 
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Fig. 1. Out-of-plane thermal conductivity k c of 
Lai. 8 _ x Euo.2Sr x Cu0 4 (x = 0, 0.08, 0.15, 0.2) as a func- 
tion of temperature. 



electronic contributions, n c is purely phononic for all Sr- 
contents. This is concluded from the very low out-of-plane 
electrical conductivity and the negligible magnetic cou- 
pling along the c-axis [5*1151110111 1| . A usual low-T phonon 
peak is present in k c for all doping levels. The gradual 
suppression of the peak upon Sr-doping is straightfor- 
wardly understood in terms of scattering by impurities 
which are induced by the Sr-ions. For x < 0.15 a step-like 
decrease is present at Tlt ~ 130 K which may be seen 
more clearly in Fig. |2 where the high-T behavior of k c is 
shown separately for each compound. At Tlt a structural 
phase transition occurs between the so-called LTO- (Low- 
Temperature-Orthorhombic, T > Tlt) and LTT-phases 
(Low- Temperature- Tetragonal, T < Tlt) - Soft tilting modes 
of the CuOg-octahedra enhance the scattering of acoustic 
phonons in the LTO-phase, which generates the step at 
Tlt- 

As revealed from Fig. [2 the T-dependence of k c in the 
LTO-phase apparently evolves from (x — 0, 0.08) an al- 
most constant behavior with a slight decrease at high T 
into a clear increase with increasing T (x = 0.15, 0.2). This 
obviously conflicts with the expected decrease of K p h (ide- 
ally following T^ 1 ). Strong deviations from the expected 
behavior could in principle arise from structural peculiar- 
ities of the LTO-phase. This can, however, be ruled out 
since in the case of x = 0.2 the LTT-, the LTO- and in ad- 
dition the so-called HTT-structures (High- Temperature- 
Tetragonal) are successively passed through with rising T 
without any influence on the increase of k c . Furthermore, 
K p h of the LTO-phase should never be equal to and never 
exceed K p h of the LTT-phase, where additional phonon- 
scattering by the aforementioned soft tilting modes does 
not exist. Yet this is the case for x — 0.15 and x = 0.2. 
For x — and x — 0.08, a similar situation is present 
concerning a hypothetical K p h of the LTT-phase which is 
extrapolated into the LTO-phase (solid lines in Fig. 0). 
These extrapolations have been obtained by fitting K p h of 
the LTT-phase with k = a/T + b, b > and extrapolating 
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Fig. 2. k c of Lai. 8 -xEuo.2Sr x Cu0 4 (x = 0, 0.08, 0.15, 0.2) 
as a function of temperature (open circles). The solid lines 
extrapolate the T-dependence of k c in the LTT-phase to high 
T (cf. text). 



this function to high T. In all cases these extrapolations 
are considerably smaller than the measured n c just slightly 
above Tlt- 1 We therefore conclude that in addition to the 
usual heat conductivity by the acoustic phonons, K aC st, a 
further transport channel for heat is present. Electronic 
and magnetic contribution to n can be excluded as al- 
ready mentioned above. The additional heat transport 
must therefore arise due to optical phonons. This is re- 
markable since in practice optical phonons are neglected 
in the analysis of phonon heat transport. 

It is commonly argued that optical phonons have dis- 
persion curves at high energy with almost no slope. In 
that case the number of excited phonons is small and 
these phonons do not propagate. Their ability to trans- 
port heat is therefore negligibly small. The situation is 
different in the case of doped La2Cu04 where it is well 
established that numerous optical phonons have a high 
velocity [T^ITH?] . Considerable transport of heat by such 
optical phonons should therefore be possible if their en- 
ergy is not too high to be excited at temperatures relevant 
here, i.e. at T < 300 K. In the following we therefore in- 
vestigate this possibility based on a simple kinetic model 
for the phononic heat transport. 

Following standard linearized Boltzmann-theory the 
thermal conductivity of a single phonon branch (labelled 
by i) in a crystal with cubic symmetry is given by [I4j 
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1 Note that an intersection between the measured and ex- 
trapolated curves also arises for extrapolation functions which 
decrease much weaker. 
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where v\., 1^, eL an d n L denote velocity, mean free path, 
energy and Bose function of a phonon. Frequently, a Debye- 
approximation with momentum independent mean free 
path is used, i.e. t>£ = v % and = l % . In this case Eq. ^ 
simplifies to 

v = ~ c ^r (2) 

with the specific heat c l v of that phonon branch. The 
usual acoustic phonon thermal conductivity K acs t is then 
obtained by summing up the contributions of the three 
acoustic phonon branches. Often the different branches of 
acoustic phonons are unified by employing an avergaged 
phonon velocity and the Debye specific heat of phonons. 
Even though such an estimate for K acs t is very crude, yet 
it is helpful to understand the T dependence of K acst in 
the low- and high-T limits: At low T the dominating scat- 
tering process is temperature independent scattering on 
crystal boundaries which leads to K acst oc cy oc T 3 . For 
T — > oo the specific heat becomes constant while the 
phonon mean free path is inversely proportional to the 
number phonons that generate umklapp processes, i.e., 
K acs t oc T _1 . These two limiting T-dependencies lead to 
the characteristic phonon-peak of K acs t at low-T which 
arises when umklapp scattering becomes important and 
over-compensates the increasing number of phonons. 

In order to obtain a description which takes optical 
phonon branches into account we take a summation over 
all relevant phonon branches (acoustic and optic), i.e., 

k = ^ = ~£>V^\ (3) 

i i 

where we extract the velocities v % from linear approxima- 
tions to the experimental (0, 0, g)-phonon dispersion re- 
lations as measured by inelastic neutron diffraction |13) . 
The specific heat c l v of each branch is calculated with the 
usual Debye model where an energy gap A l opt has to be 
considered for the optical branches. Hence we have 

4 = 3A^s(Jr) i\ (4) 

where Na and ks are Avogadro's and Boltzmann's con- 
stants, l D — v l hkg (G^n) 1 ' 3 with n as the density of 
unit cells and 

Ii= [^( x+ ^A 2 * 2 e X p(, + ») ^ (5) 
Jo \ T J [expOr + - l] 2 

The different v % , A^ pt and l D of the phonon branches 
which we used for our analysis are summarized in Tabled 

It is evident from the table that the v l of the optical 
and acoustic branches are of comparable magnitude and 
therefore a significant propagation of the optical modes 
is present indeed. The degree of excitation of the optical 
modes is reflected by the specific heat, which is shown in 
Fig. |3| (a), for the different phonon branches. Remarkably, 
at T = 300 K the c\ of the optical branches labelled A 



Type 


^ Pt (K) 


v % (ms" 1 ) 


Ob (K) 


acoustic 1+2 





5423 


280 


acoustic 3 





12289 


635 


optic A 


251 


4263 


220 


optic B 


333 


3108 


161 


optic C 


359 


28041 


1449 



Table 1. Optical gaps A l opt , velocities v* and 'Debye' temper- 
atures & % D of the considered phonons. 



and B already amount to about 80% of the Dulong-Petit 
value, i.e., the major part of these optical modes is already 
populated at this T? c v of the strongly dispersing branch 
C apparently evolves at much higher T. 

In our kinetic model the products c v v l , which measure 
the excitation of propagating phonons, are more relevant 
for the heat transport. In Fig. [31(b) the T-dependence of 
the sum Stotai = J2i c v yl °f all contributions is shown 
as a solid line while the different contributions 5 acs t — 
J2T St c v yt an d ^opt = X^ Pt c v v% °f acous ti c and optical 
phonons are displayed as broken lines. Apparently, prop- 
agating optical phonons become significantly excited at 
T > 70 K and already contribute 40% of <S tota ] at 300 K. 

It is now very interesting to use this information and 
analyze experimental phononic thermal conductivities in 
terms of acoustic and optic contributions, K acs t and K Q pt- 
If we follow our model, an accurate separation of K acs t and 
K op t is now possible with the separate phonon mean free 
paths l l of each branch. These are, however, not known. 
As a first approximation we therefore assume an integra- 
tive mean free path l l = ^ p h- The total phonon thermal 
conductivity then separates into acoustic and optic con- 
tributions in the same ratio as iS acs t and 5 op t . As an exam- 
ple, we have calculated K acst and K op t of Lai.sEuo.2Cu04 
as a function of T as shown in Fig. [21(c). 3 Apart from the 
remarkable high /t pt> the T-dependence of K acs t at high 
T is now indeed in accordance with the expected decrease 
with increasing T. In fact, at T > 200 K, K acs t even fol- 
lows a T~ ^dependence as is illustrated in Fig. [3] (c) by 
the solid line which represents a fit K acst = a/T + b, b > 0. 

The inset of Fig. |3] depicts the T-dependence of / p h of 
Lai.gEuo.2Cu04 which we have calculated from our ex- 
perimental data and 5tot a i- Its strong decrease has to be 
attributed to the growing importance of umklapp scatter- 
ing processes with increasing T. Note that Z p h ~ 15 A at 
T w 350 K. This value is slightly larger than the lattice 
constant c = 13.2 A which defines the lower limit for i p h. 
At high T, our estimation of K opt is therefore close to a 
lower limit of K op t, i-e., optical contributions could even 
be more important than it is revealed in our study. 

Several simplifications in our approach should be men- 
tioned besides the integrative l p h- First, the linear ap- 
proximation of the dispersion curves neglects curvatures 

2 Since only a single phonon branch is considered, the 
Dulong-Petit value is 24.9/3 Jm^K" 1 = 8.3 Jm -1 K -1 . 

3 As a first approximation we assume for the calculation of 
/tacst and K opt that the structural phase transition at Tlt ~ 
130 K does not affect the lattice dynamics. 



4 



Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 




Temperature (K) 

Fig. 3. (a) Debye-specific heats c v of acoustic and optic 
phonon branches as a function of T. (b) Different summations 
(optical, acoustic and acoustic+optic) of the products v l c v as 
a function of T. (c) k c (full circles) of Lai.sEuo^CuCU and the 
calculated contributions by acoustic acoustic (open squares) 
and optic (open triangles) contributions. The solid line illus- 
trates the T _1 -dependence of K acst at high T. Inset: l p h of 
Lai.8Euo.2CuC>4. 



of the real dispersion relations which are present in some 
parts of the Brillouin-zone. Second, even though in the 
real material comparable dispersing optical phonon modes 
are present along all crystal directions, it is certainly not 
perfectly isotropic as implied by our approach. Third and 
finally, apart from the three acoustic branches we only 
take three optical branches into account. It is not clear 
whether the first two simplifications influence our estima- 
tion towards a smaller or larger result for K pt- Yet it is 
unlikely that the order of magnitude of At opt is significantly 



smaller. Taking more optical branches into account, as re- 
quired for a more realistic description, will certainly shift 
K op t towards higher values. 

Our data suggest a clear doping dependence of the 
ratio of K op t and K aC st since the high-T increase of k c be- 
comes more and more apparent with increasing Sr-content 
(cf. Fig. |2J) • We note in this context that for some of the 
dispersive optical phonon branches the optical gaps and 
the slopes of the dispersion curves change upon Sr-doping. 
Therefore, a growing importance of K op t with increasing 
Sr-content could, for example, be related to decreasing 
optical gaps and increasing slopes of the dispersion curves 
[12113) . However, in order to verify such a detailed inter- 
pretation, a more accurate description of K p h than our 
simple approach is necessary. In principle, such can be 
achieved when the full set of «£, and l^ = v^t^ (with 
the phonon life time t£), accessible through neutron scat- 
tering experiments, is employed into the >r 4 in Eq. ^ 

In conclusion, we have used a simple model in order 
to analyze the phononic out-of-plane thermal conductivity 
k c of Lai.g_ x Euo.2Sr x Cu04 in terms of contributions by 
acoustic and optic phonons. Heat conductivity by optical 
phonons turns out to be significant in this material, since 
according to our approach it amounts at room tempera- 
ture to about 40% of the total phonon thermal conduc- 
tivity. This result has important implications for studies 
on magnetic heat transport of low dimensional spin sys- 
tems which are currently on debate [Tj;2,3.4,5,15.16|l7|. 
Unusual T-dependencies of k which are hastily attributed 
to magnetic heat transport must carefully be excluded to 
stem from heat transport by dispersive optical phonons. 
We stress that more careful analysis of this usually under- 
estimated mechanism of heat transport is desirable. 

C.H. acknowledges support by the DFG through grant 
HE3439/3-1. We thank M. Braden for helpful discussions 
and N. Jenkins for proofreading the manuscript. 
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